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Abstract 

We study the structure and non-perturbative properties of a special Green's function, u{q 2 ), 
whose infrared behavior has traditionally served as the standard criterion for the realization of 
the Kugo-Ojima confinement mechanism. It turns out that, in the Landau gauge, u(q 2 ) can be 
determined from a dynamical equation, whose main ingredients are the gluon propagator and 
the ghost dressing function, integrated over all physical momenta. Using as input for these two 
(infrared finite) quantities recent lattice data, we obtain an indirect determination of u{q 2 ). The 
results of this mixed procedure are in excellent agreement with those found previously on the 
lattice, through a direct simulation of this function. Most importantly, in the deep infrared the 
function deviates considerably from the value associated with the realization of the aforementioned 
confinement scenario. In addition, the dependence of u(q 2 ), and especially of its value at the 
origin, on the renormalization point is clearly established. Some of the possible implications of 
these results are briefly discussed. 

PACS numbers: 11.15.Tk 12.38.Lg, 12.38.Aw, 
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I. INTRODUCTION 



The problem of quark confinement and gluon screening is of central importance in QCD, 
and a large body of work has been dedicated to its understanding [1]. Some of the most 
widely explored mechanisms attempting to explain how quarks confine make concrete predic- 
tions about the non-perturbative behavior of the fundamental Green's functions of the the- 
ory. For example, a central ingredient in the center vortex picture of confinement put forth by 
Cornwall is the dynamical generation of a gluon mass [2] through the well-known Schwinger 



mechanism jj, 4,|5j, implemented within the pinch technique (PT) framework 
addition to taming the infrared divergences intrinsic toperturbation theory (Landau pole), 
this mass gives rise to a low energy effective theory [8j which supports quantum solitons 
(center vortices), not present in the massless theory, whose condensation furnishes an area 
law to the fundamental representation Wilson loop, thus confining quarks On the other 
hand, the adjoint potential shows a roughly linear regime followed by string brea king when 



aaa. 



In 



the potential energy is about 2m, where m is the induced mass of the gluon [lOl . Illl ] corre- 



sponding to gluon screening 



12 



13( . At the level of the two fundamental Green's functions 



of the theory, namely the gluon and ghost propagators, the predictions of the above picture 
are very definite: the gluon propagator is infrared finite (due to the generation of the gluon 



mass 



15], whose phenomenological value has been delimited in Ref. 14j), while, as has been 



shown recently 



15 



16), in the Landau gauge the ghost remains massless, but with a finite 



dressing function (due to the saturation produced by the gluon mass) [171 ]. 

An entirely different set of predictions is obtained within the Kugo-Ojima (KO) scenario, 
which also establishes a highly non-trivial link between confinement and the infrared be- 
havior of some of the most fundamental Green's functions of QCD. In the KO confinement 
picture (in covariant gauges), the absence of colored asymptotic states from the physical 
spectrum of the theory is due to the so-called quartet mechanism [18]. A sufficient condition 
for the realization of this mechanism (and the meaningful definition of a conserved BRST 
charge) is that a certain correlation function, to be denoted by u(q 2 ), defined in Eg. (12.211) . 

2oj, in the 



19|. In addition, as first noted by Kugo 



should satisfy the condition u(0) = - 
Landau gauge, w(0) is related to the infrared behavior of the ghost dressing function F(q 2 ) 
[see Eq. (12. 7p ] through the identity F _1 (0) = 1 + u(0). Therefore, the KO confinement sce- 
nario predicts a divergent ghost dressing function, and vice-versa. Interestingly enough, the 
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same prediction a 
horizon condition 



)out F 1 (0) is obtained when implementing the Gribov-Zwanziger (GZ) 



21 



221 ] : in the IR region the ghost propagator diverges more rapidly than 



at tree-level [231 ]. Furthermore, it has been also argued that the Landau gauge gluon prop- 



agator should vanish in the same limit 



24|. This alleged connection between confinement, 



the horizon condition, and an infrared "enhanced" ghost dressing function has served as the 
theoretical cornerstone of the so called "ghost- dominance" picture of QCD 251 ] . 



Turns out that recent large volume lattice simulations [26j, [27[ appear to be at odds 
with the original KO and GZ pictures described above, at least as far as their predictions 
about the infrared behavior of the Green's functions are concerned 2a] . Specifically, various 
lattice studies, both in SU (2) and SU (3), find (in the Landau gauge) an infrared finite gluon 



propagator 



29| | and an infrared finite ( "non-enhanced" ) ghost dressing function. Evidently, if 



taken at face value 



301 ] , these results furnish strong support for the PT picture of dynamical 



gluon mass generation and the ensuing confinement mechanism. 

It is perfectly clear that further detailed scrutiny from all possible angles must be imple- 
mented before reaching a definite conclusion on any of these issues. In this vein, it is natural 
to ask what one really knows about the KO function u(q 2 ). Turns out that u(q 2 ) has been 



studied directly on the lattice by Sternbeck 



31], using the field-theoretic definition of u(q 2 ) 



appearing in the KO formulation. As can be plainly seen from Sternbeck's result (repro- 
duced for convenience in Fig. [8] of this article), u(0) deviates appreciably from its KO value 
of —1; specifically, the function u(q 2 ) saturates in the deep infrared around approximately 
—0.6. Interestingly enough, in a recent article 32j Kondo gave a simple derivation of this 
same value, after appropriately modifying the KO construction in order to self-consistently 
accommodate the GZ horizon condition. 

Quite remarkably, in the (background) Landau gauge [33J] the KO function coincides 
with a certain auxiliary function, usually denoted by G(q 2 ), which constitutes a crucial in- 
gredient in the modern formulation of the PT by means of the Batalin-Vilkovisky (BV) 
quantization formalism 34j. Specifically, G(q 2 ) is the form-factor multiplying g^ v in the 
Lorentz decomposition of a special Green's function, denoted by A Mi ,(g), which enters in 



all "background-quantum" identities 



35 



361 ]. i.e., the infinite tower of non-trivial relations 



connecting the Green's functions of the background field method (BFM) 37| to the con- 
ventional ones (e.g. gauges). Notice also that G(q 2 ) plays a prominent role in the new 
Schwinger-Dyson (SD) equations derived within the PT framework |38f| , which, due to the 
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special properties of the Green's functions involved, can be truncated in a manifestly gauge 



invariant way 



39|. 



As has been shown in a recent article 



401 ]. one may derive a dynamical (SD-like) equation 



for G(q 2 ), which, under mild assumptions, allows one to reconstruct G(q 2 ) from the knowl- 
edge of the gluon and ghost propagators. Specifically, G(q 2 ) is determined by integrating 
over all virtual momenta (k) a kernel involving the product A(k)F(k + q). We empha- 
size that the aforementioned dynamical equation is not a simple relation of several Green's 
functions at some special isolated point; instead, the value obtained for G(q 2 ) at any point 
(such as q 2 = 0) must be compatible with the behavior of F and A in the entire range of 
their physical (euclidean) momenta. In particular, one must know their behavior not only in 
the IR, but also in the intermediate region of momenta (0.3-3 GeV), which appears rather 
difficult to obtain from SD studies [411 ] . This feature is very powerful, because it probes the 
details of the fundamental Green's functions over an extended range of momenta, rather 
then just a single point. 

In the present work, we use the available lattice data on the gluon and ghost propagator 
as input into the aforementioned dynamical equation, thus obtaining an indirect determi- 
nation of G(q 2 ) in the entire range of available lattice momenta. Given the Landau gauge 
coincidence between G(q 2 ) and u(q 2 ), this procedure automatically determines the KO func- 
tion as well. This, in turn, permits us to obtain the value of the KO parameter w(0), as well 
as the GZ horizon function, and study their dependence on the renormalization point fi. 
Our analysis reveals an impressive self-consistency between the various ingredients entering 
into the calculation. In particular, the results obtained through our combined method (SD 
using lattice data as input) are in excellent agreement with those of 3J|, obtained through 
a direct simulation of the KO function. 

The paper is organized as follows. In Section [II] we briefly introduce the BV framework 
for SU(N) Yang-Mills theories, where the function A M! ,(g) appears naturally Next, we 
review a number of relations where this function plays a key role: (i) the background- 
quantum identity relating the conventional and the BFM gluon propagators; (ii) the relation 
between the ghost dressing function F(q 2 ) and the A M „(g) form factors G(q 2 ) and L(q 2 ); 
(in) we establish the crucial equality u(q 2 ) = G(q 2 ); (iv) the relation with the GZ horizon 
function. The central results of this article are presented in Section IIHI Specifically, the 
Lorentz decomposition of A Mi ,(g) gives rise to two form-factors, the G(q 2 ), which in the 
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previous section has been identified with the KO function u(q 2 ), and the L(q 2 ), which has the 
particular property of vanishing in the deep IR. After establishing the dynamical equations 
governing G(q 2 ) and L(q 2 ), we use the recent lattice data for the gluon and ghost propagators 
as input in these equations. In addition to the equations for G(q 2 ) and L(q 2 ), we consider 
the SD equation for the ghost, which is calibrated in order to be numerically compatible 
with the lattice data (at an impressive precision) simply by adjusting the gauge coupling 
to values that are slightly above the standard two-loop MOM prediction; the obtained 
value of the coupling is then used into the equations for G(q 2 ) and L(q 2 ). We use the 
multiplicative renormalizability of the gluon and ghost propagators in order to rescale the 
lattice data to different values of the renormalization point. Even though this procedure has 
an intrinsic limitation set by the relatively short reach of the available data into the UV, it 
amply demonstrates that u(q 2 ) depends non-trivially on /i, in excellent agreement with the 



observation established in 



31] . Finally, in Section [TV] we present our conclusions. 



II. CONNECTING THE KUGO-OJIMA AND G FUNCTIONS 

As already mentioned in the Introduction, in the Landau gauge the KO function u(q 2 ) 
may be shown to be identical to the function G(q 2 ), which appears in several formal contexts. 
In this section we first formulate Yang-Mills theories in the BV framework, which allows the 
derivation of a tower of identities, whose common ingredient is the function G(q 2 ). Then 
we will show why u(q 2 ) = G(q 2 ), and will review the connection between u(q 2 ) and the GZ 
horizon. This main purpose of this section is to serve as a reminder and to bring together 
various seemingly disjoint pieces of information. For specific details on each topic the reader 
is referred to the corresponding extensive literature. 



A. Batalin-Vilkoviski formalism 



In the BV formulation of Yang-Mills theories 34j, one starts by introducing certain 
sources (called anti-fields in what follows) that describe the renormalization of composite 
operators; the latter class of operator is in fact bound to appear in such theories due to the 
non-linearity of the BRST transformation of the elementary fields. In much the same way, 
the quantization of the theory in a background field type of gauge requires, in addition to 
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the aforementioned anti-fields, the introduction of new sources which couple to the BRST 



variation of the background fields 



35] . These sources are sufficient for implementing the full 



set of symmetries of a non-Abelian theory at the quantum level, and in the case of quarkless 
SU(N) QCD, lead to the master equation 



5T 5T 5T 5T om ST om / 5T 5T 
+ - — -r— + B m -— + Q m 



5 A*" 1 bA™ Sc m 5c m Sc m M SA m 5 A™ 



0. (2.1) 



In the formula above, T is the effective action, A* and c* the gluon and ghost anti-fields, A is 
the gluon background field, and Q the corresponding background source; finally B denotes 
the Nakanishi-Lautrup multiplier for the gauge fixing condition. 

To determine the complete algebraic structure of the theory we need two additional 
equations. The first one is the Faddeev-Popov equation, that controls the result of the 
contraction of an anti-field leg with the corresponding momenta. In position space, it reads 

'" + ^TI7 -(2>*n M r = 0, (2.2) 



6c m \ 5A* 



where {V^) m = <9^$ m + gf mnr A n ^ r [in the case of (2^$) m replace the gluon field A with 
a background gluon field A]. The second one is the anti-ghost equation formulated in the 
background field Landau gauge, which reads [33 ] 

( V» I - (V^A*) m - f mnr c * n c r + f mnr c r = 0, (2.3) 

5c m V 5B n v ; 

This equation fully constrains the dynamics of the ghost field c, and implies that the latter 
will not get an independent renormalization constant. The local form of the anti-ghost equa- 
tion (I2.3P is only valid when choosing the background Landau gauge condition (V^A^)" 1 = 0; 
in the usual Landau gauge, d^A™ = 0, an integrated version of this equation is available. In 
fact, even though the results that follow will be derived for convenience in the background 
Landau gauge, they are valid also in the conventional Landau gauge of the R^. 

Now, differentiation of the functional fl2.ll) with respect to a combination of fields con- 
taining at least one ghost field or two ghost fields and one anti-field (and setting the fields 
and sources to zero afterwards) will provide the Slavnov- Taylor identities of the theory. 
Differentiation with respect to a background source and background or quantum fields will 
provide, instead, the so called background-quantum identities [35j, |36[], which relate 1PI 
Green's functions involving background fields with those involving quantum fields. Finally, 
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k, a 




H m (k,q) = if™ + „ - ■ 

q — k 



FIG. 1: Diagrammatic representation of the functions A and H. 



differentiation of (12 .2p and (12.31) with respect to fields and anti-fields or background sources 
give rise to relation among the different auxiliary ghost functions appearing in the theory. 

The important point is that, when carrying out these differentiations, the following func- 
tion appears (Fig. [Q 

iA^(g) = r nM A*(?) 

= g 2 C A [ H^D(k + q)A^(k)H (n/ (k,q), 
Jk 



and (in <i-dimensions) 
G(q 2 ) 1 



^G(g 2 )+^L(g 2 ), 



1 



(2.4) 



{dq^A, u - g 2 A|;) . (2.5) 



(rf-i)g 2 ^ " * * v " y (rf-i)g 2 

In the equations above, the color factor S mn has been factored out (as always in what follows), 
Ca represents the Casimir eigenvalue of the adjoint representation [Ca = N for SU(N)), 
and j k = // 2£ (27r) J d d k, with d = 4 — e the dimension of space-time. and .D represents 
the gluon and ghost propagator respectively, defined as 



2^ . . iF(Q 



P^(g)A(g 2 )+£ 



D(q 2 



(2.6) 
(2.7) 



where £ denotes the gauge-fixing parameter, and P^ u (q) = Q\n> — Q^Qu/q 2 is the usual trans- 
verse projector; notice that A _1 (g 2 ) = q 2 + ill(q 2 ), with H^ v (q) = P^ u (q)Tl(q 2 ) the gluon 
self-energy. F(q 2 ) is the so called ghost dressing function. Finally, the function H^ v {k,q) 



(see Fig. [T] again) is in fact a familiar object, since it appears in the all-order Slavnov- Taylor 



identity satisfied by the standard three-gluon vertex 



42] 



q a T afiV {q,k h k 2 ) = iF^A-^kDP^H^h) 

- iF^A-^kDP^H^h). (2.8) 

It is also related to the full gluon-ghost vertex T^(k, q) by the identity 

q v H^(k,q) = -iT„(k,q). (2.9) 

At tree-level, Hfj = ig^ u and T^\k,q) = T^(k,q) = —q^. 

B. Background-quantum identities 

The first identities where the function A M!/ appears are the so-called background-quantum 
identities, i.e., the infinite tower of non-trivial relations connecting the BFM Green's func- 
tions to the conventional ones 



35 



36( | . Consider, for example, the result of differentiating 
the functional ( 12. II) with respect to a background source and a background gluon, on the 
one hand, and a background source and a quantum gluon, on the other, one obtains two 
equations, 

Using the transversality of the gluon two-point function, these two equations can then be 
appropriately combined to yield the important identity 

r^ Av (q)=[l + G(q 2 )} 2 T A ^(q), (2.11) 

or, in terms of propagators 

A- 1 (q 2 )=[l + G(q 2 )} 2 A- 1 ( q 2 ). (2.12) 

The quantity A(g 2 ) appearing on the left-hand side of the above equation captures the 
running of the QCD f3 function, exactly as happens with the QED vacuum polarization; 
this is a fundamental property of the BFM gluon self-energy, valid for every value of the 
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(quantum) gauge-fixing parameter [37|. This can be easily checked to lowest order, where 
Eqs. (12. 4p and (12. 5p give (in the Landau gauge) 



1 + GV) = 1 + 7 



4 48tt 2 



and thus 



9 



2 48tt 2 



2\1 



(2.13) 



(2.14) 



where 6 = 11Ca/487t 2 is the first coefficient in the QCD f3 function. Eq. (12.121) plays a 
central role in the derivation of a new set of SDEs [38] that can be truncated in a manifestly 
gauge invariant way [391 ] . 



Let us conclude this subsection by noticing that in more general identities the function 
L(q 2 ) is also relevant. Consider e.g., the identity relating the three-gluon proper vertices. 
One then has 



r A tl A a A ( k ^ fc 2) = [gfc + A£(g)] r AvAaAp (k u k 2 ) + 



[1 + G(q 2 )} T A ^ AaAp (k u k 2 ) + q J^L{q 2 )T AvAaAf) {k u k 2 ) + ■ ■ ■ , (2.15) 



where the omitted terms involve other auxiliary Green's functions (see 35|, |36j), irrelevant 
to our discussion. 



C. Two-point ghost sector 



Let us now consider the two-point functions. Differentiating the ghost equation (12. 2p 
with respect to a ghost field and a background source we get the relations 



Fcc(q) = -iq u T cAt (q), 
T cn^(q) = q fi + q u A ll u(q)- 



(2.16) 



On the other hand, differentiating the anti-ghost equation (12.31) with respect to a gluon 
anti- field and an anti-ghost, one gets 



T cA t (q) 

r c5 (g) 



q u + q f *A flv (q), 



(2.17) 



/i J -V /i u w 



FIG. 2: Connected components contributing to the function G™™(q). 



Next, contracting the first equation in (I2.17P with q u , and making use of the first equation 
in (12.161) . we see that the dynamics of the ghost sector is entirely captured by the A^ 
auxiliary function, since 

iT c ^q) = q 2 + q^q u A^{q). (2.18) 
Introducing the Lorentz decompositions 



r c A*(g) = qf,C(q 2 ), Tg^O) = q tl E{q 7 



(2.19) 



we find that Eq. (12.181) together with the last equation of (12.161) and (12. 17ft give the identi- 



ties 



20 



33] 



C(q 2 ) = E(q 2 ) = F-\q 2 ), 
F-\q 2 ) = l + G {q 2 ) + L{q 2 ). 



(2.20) 



Finally, recalling that the dimension of the gluon anti-field A* is three, while the dimension 
of the Q source is one, power counting shows that (i) all functions appearing in Eqs. (I2.16P 
and (12.171) a re divergent, and (ii) the divergent part of A^ u (q) can be proportional to g^ u 



only 



33. 



40J 



D. The (background) Landau gauge equality between u(q 2 ) and G(q 2 ) 

A crucial ingredient to our analysis is the equality between the KO function and the 
G(q 2 ), in the (background) Landau gauge [33}. To see this, we start with the following 
(Euclidean) two-point function of composite operators 

J d 4 x e-^-^(T[(2V)™ {f nrs A n v c s ) y ]) = P^(q)5 mn u(q 2 ), (2-21) 

which, due to the identity 

J d 4 x e- iq - (x ~ y) (T [ (2V)™ c n ] ) = -i^5 mn (2.22) 

10 



can be related to the function ([ (U^c)™ (X^c)™] ), through 

J d 4 x e-^-^(T[(2V)™ (V.c);}) = -M^™ + P^(q)5 mn u(q 2 ). (2-23) 

On the other hand, observe that in the background Landau gauge the function appearing 
on the lhs of the above equation is precisely given by 

- e " (9 > = «f^' (2 - 24) 

where W is the generator of the connected Green's functions, and the two connected dia- 
grams contributing to Q^ v are shown in Fig. [2j Factoring out the color structure and making 
use of the identities (12.201) one has 

-iGM = ^u(q) + T n ,-c(q)D(q 2 )T Atc (q) 

= -^r + P»M)G(q 2 ). (2.25) 

Passing to the Euclidean formulation, and comparing with Eq. (I2.23f) . we then arrive at the 
important equality 

u{q 2 )=G{q 2 ). (2.26) 

Then, the usual KO confinement criterion may be equivalently cast in the form: 1+G(0) = 0. 
Evidently, if L(0) = [see discussion after Eq. (13.61) ]. then from the identity (I2.20p follows 
that if the KO criterion is satisfied then the ghost dressing function diverges in the IR. 

E. Gribov-Zwanziger horizon 



In order to avoid Gribov copies 2l|, in the GZ formulation of Yang-Mills theories the 



partition function assumes the form (in (^-dimensional Euclidean space) [2 
z i = J [dA]5(d fM A fl )DetM exp |-S YM + 1 J d d x h(x) 



(2.27) 



where Sym is the Yang-Mills action, M = —d^V 1 is the Faddeev-Popov operator, and the 
functional h(x) = h[A](x) is the so-called GZ horizon function given by 

h(x) = - jd d ygr mr A^(x)(M- 1 Y s (x,y)gr n A:(y) ] (2.28) 
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thus, the action corresponding to the partition function above clearly contains a non-local 
term. The GZ parameter, 7, is determined through the so-called horizon condition, which 
for SU(N) assumes the form 

(h( x ))~/ — d(N 2 — 1). (2.29) 

This condition can be rewritten in terms of the vev of the GZ horizon function if we integrate 
both sides over d d x. In this way we get 



= L [d d x (h(x))„ = d(N 2 - 1). (2.30) 

Vd J 

On the other hand, assuming that 7 is small, one can expand in powers of 7; retaining the 
first order only, one gets 

(MO)) 7 ^(MO)) 7=o + 0(7)- (2.31) 

The right-hand side of the above equation can be related to the trace of the following 
Green's function (Euclidean space) 



H™(q) = 6™ 

n 

in the limit q — >• [32]. Specifically, 



p p M + 9j f F (Q 2 ) 



KM, (2-32) 



(M0)) 7=0 = ^ [d d x(h(x)) 7=0 = -limTr{^(g)} 

= -(N 2 - 1) {(d - 1)G(0) + F(0) [G(0) + L(0)]} . (2.33) 

This result allows to rewrite the GZ horizon condition in terms of G(0) (and therefore 
of the KO parameter w(0)); this will, in turn, restrict the allowed values of u(0). In the 
limit of vanishing Gribov parameter, one can use the result of f !2.33|) to solve the horizon 
condition (I2.29p . in the approximation (j2.3ip ; If -^(0) = 0, one finds (in d = 4) the following 



value of the KO parameter 



32] 



w(0) = G(0) = -| (2.34) 



which is very close to that obtained directly from the lattice [3jJ , and, as we will see in the 
next section, from our independent analysis. 
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III. EXTRACTING THE KUGO-OJIMA FUNCTION FROM THE LATTICE 



In this section we study the behavior of the function G{q 2 ) by using the available lattice 
data on the gluon and ghost propagators. Specifically, we will first write down the dynamical 
equations that govern the functions G(q 2 ) and L(q 2 ), which involve both A(q 2 ) and D(q 2 ), 
and using as an input the lattice results for these propagators we will obtain an indirect 
determination of G(q 2 ). Of course, by virtue of the fundamental equality u(q 2 ) = G(q 2 ), 
any information on G(q 2 ) translates automatically to the Kugo-Ojima function. 



A. Dynamical equations 

The dynamical equations governing G(q 2 ) and L(q 2 ) may be obtained directly from the 
defining equation (12.4j) . by appropriately contracting it and taking its trace. Specifically, 
from Eq. (12.41) one has (in <i-dimensions) {40} 



G { q 2 ) = g± 
L(q 2 ) - ! '~ C> 



d 



[ AP°(k) H ap (k, q)D(k + q) + i\ [ q p A pa (k) T°(k, q)D(k + q) 
Jk <? Jk 

4 / q P A pa (k) T°(k, q)D(k + q)+f A^(k) H ap (k, q)D(k + q) 
. Q Jk Jk 



.(3.1) 



Adding the above equations by parts, and employing (12.91) . one may easily demonstrate |40j], 
in the Landau gauge, the validity of the second identity in Eq. (12.201) . given that the standard 
SD equation for the ghost propagator (Fig. [3]) reads 

iD-\q 2 ) = q 2 + ig 2 C K [ T p A pu (k)T"(k, q)D(q + k). (3.2) 

Jk 

The two equations in (13.11) involve five basic ingredients: the two-point functions A(q 2 ) 
and D(q 2 ), the vertex functions r p (k,q) and H^ v (q, k), and, eventually, the value of the 
(renormalized) coupling g 2 , at different renormalization points. Knowledge (direct or indi- 
rect) of these ingredients (for example from the lattice) would, in turn, determine fully the 
functions G(q 2 ) and L(q 2 ). 

As far as the vertex functions T p (k, q) and H^k, q) are concerned, due to the transver- 
sality of the gluon propagator, the only relevant form factors are 

-T p (k,q) = B 1 (k,q)q tl , 
-iH^(k,q) = A 1 (k,q)g pu + A 2 (k,q)q^q u . (3.3) 
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( > O > ) 1 = 

q 



+ 




k + q 



FIG. 3: The SDE satisfied by the ghost propagator. 

In order to study the relevant equations further, we will approximate the form factors above 
with their tree-level values, i.e., Ai(k,q) = B\(k,q) = 1, and A 2 (k, q) = 0; according to 



lattice studies 
(13. ip become 



431 ]. this appears to be a very good approximation. Then, the equations in 



G(q 2 ) 
L{q 2 ) 



9 2 C A 
d- 1 
9 2 C A 
d — 1 .i k 



(d - 2) + 



(fc.g) 2 



A(A;) J D(fc + g), 



1 - d 



(fc • g) ; 

k 2 q 2 



A(k)D(k + q). 



(3.4) 



Now, as discussed in detail in 



40] the (unrenormalized) equations in (13.41) must be prop- 



erly renormalized, i.e., in such a way as to preserve the validity of the (BRST-induced) 
second identity in ( 12.201) . This finally amounts to renormalizing these equations through 
multiplication by Z c , the renormalization constant of the ghost propagator. Specifically, 
after introducing 

A(gV) = Z A V)A (g 2 ), 
F(gV) = Z-\^)F (q 2 ), 



#0 2 ) = Z \fi 2 )g , 



we obtain in the Euclidean space, after setting q 2 
implementing the standard angular approximation, 



x, k 2 



(3.5) 



y and a s = g 2 /4:7i, and 



1 + G(x) 
L(x) 



12tt 



167T 

F(x) 



^ fdy V (3 + f) A(y) + [°°dy U + f) A( y )F(y) 
x Jq \ ox/ j x \ 6y J 

Jx y 



(3.6) 



From this last equation it is easy to see (e.g., by means of the change of variables y = zx) 
that if A and F are IR finite, then L(0) = 0, as mentioned before 441 ] . 

After imposing the renormalization condition F(fi 2 ) = 1, one finds that the renormaliza- 
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Gluon Propagator 

L=64, p=5.7 
» n = 3.0 GeV 
• (i = 3.6 GeV 
H = 4.3 GeV 

fit n = 3.0GeV 

fit n = 3.6 GeV 

fit n - 4.3 GeV 




q [GeV] 

FIG. 4: Lattice results for the gluon propagator renormalized at three different renormaliza- 
tion points: fj, = 3.0 GeV (black curve), fi = 3.6 GeV (red curve) and [i = 4.3 GeV (green 
curve). We also show the corresponding fits using Eq. (|3.1U|) . The fitting parameters are: 
a = 0.162GeV 2 , b = 0.367 GeV -1 and c = 1.5 (/x = 3.0GeV); a = 0.147GeV 2 , 6 = 0.334GeV" 1 
and c = 1.5 (/u = 3.6 GeV); a = 0.137 GeV 2 , b = 0.311 GeV" 1 and c = 1.5 (/i = 4.3 GeV). 

tion constant Z c appearing in the equation for G(x) is given by 
a s C A 



1 



16tt 



(3.7) 



This renormalization constant renormalizes Eq. (13.21) as well, when the latter is evaluated 
under the same approximations, i.e., 

a s C A 



X 



Z r - 



167T 



x Jo v xJ 



dy ( 3 - - ) A(y)F(y) 



(3.8) 



Note that one cannot choose simultaneously the condition G(fi 2 ) = once F(fi 2 ) = 1 has 
been imposed; indeed, given that L(^i 2 ) ^ 0, such a choice would violate the identity of 
Eq. fT2T20|) . 

B. Numerical analysis 



The starting point for our numerical analysis are the lattice results for the gluon prop- 
agator A(g 2 ) reported in 27]. In order to eventually study the dependence of the KO 
function on the renormalization point, we would like to obtain the lattice data at different 
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0,5-| , , , 
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q 2 [GeV 2 ] 




FIG. 5: Top left panel: The ghost dressing function F(q 2 ) obtained from SDE (black continuous 
line) compared to the lattice data of [stJ at (j, = 3.0 GeV. Top right panel: Same as in the previous 
panel but renormalized at \x = 3.6 GeV. Bottom left panel: Same as before but renormalized at 
fx = 4.3 GeV. Bottom right panel: The SDE solutions for the three different renormalization points 
all together. 

renormalization points. Since the gluon propagator is multiplicatively renormalizable, the 
relation [4q] 

a, 2 2, A(gV 2 ) 

can be used to connect a set of points renormalized at n with the corresponding set renor- 
malized at v. Choosing the three different values fi = {3.0, 3.6, 4.3} GeV, we then obtain 
the three curves shown in Fig. HI In the range of available momenta, a very accurate fit is 
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3.0 GeV and a(n ) = 0.388 
3.6 GeV and atf) = 0.330 
4.3 GeV and <i( M 2 ) ; 0.295 




1E-3 0,01 



q[GeVl 



"l I 

- |i = 3.0 GeV and afci 2 ) = 0.388 

- ix = 3.6 GeV and ctfc 2 ) = 0.330 
n ; 4.3 GeV and g(n') ; 0.295 




q [GeV] 



FIG. 6: Left panel: — G{q 2 ) determined from Eq. (13.6|) . using the solutions for A(q 2 ) and D(q 2 ) 
presented in Figs. H] and [5] at the same renormalization points [i. Right panel: Same as in the 
previous panel but this time for L{q 2 ). 



provided by the expression 



A(g 2 



(3.10) 



a + bq 2c ' 

as shown by the continuous line in Fig. H] (the values of the fitting parameters a, b, and c 
are also reported there). 

The next step is to employ the ghost SDE given in (13.81) in order to deduce the appropriate 
values that one must use for the gauge coupling. To that end we will follow the following 
steps: (i) employing once again the relation (13.91) . with A — > D, we generate from the lattice 
data on F(q 2 ) reported in 27j the data sets for F(q 2 ) corresponding to the renormalization 
points used previously for A; (ii) using as input in (13. 8ft the different sets of results obtained 
in the previous step for A, we solve the integral equation (13.81) numerically, thus determining 
F(g 2 ;/i 2 ); (hi) the values of a s (/i 2 ) are fixed by demanding that the solutions obtained in 
(ii) match the different lattice sets generated at step (i). 

The results of this procedure are displayed in Fig. [5j where we show both the com- 
parison between the lattice data and the solutions of Eq. (13.81) (first three panels), as 
well as the dependence on /z 2 of these solutions (fourth panel). The couplings found are 
a(fi 2 ) = {0.388, 0.330, 0.295} for p = {3.0, 3.6, 4.3} GeV respectively. Note that the val- 
ues of a; s (/i 2 ) obtained by this procedure are about 20% higher than those found from the 



two-loop MOM calculation of 



46). 
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RGI product 

a(n 2 ) = 0.388 and n = 3.0 GeV 

— = 0.330 and i± = 3.6 GeV 
a(n 2 ) = 0.295 and n = 4.3 GeV 




q[GeVl 



FIG. 7: The renormalization-group invariant product d(q 2 ) obtained combining our results for 
A(q 2 ) and G(q 2 ) according to Eq. (pTTTj) . 

At this point, all necessary ingredient for determining the functions G(q 2 ) and L(q 2 ) are 
available. Substituting them into the corresponding equations given in (13.61) . we obtain the 
solutions shown in Fig. [6j Notice that L(q 2 ) vanishes in the deep IR, as expected. 

A very stringent test of the quality of the obtained solutions can be devised by observing 
that, on formal grounds, the combination 



2 <rV)A(g 2 ;/^ 



[1 + G(q 2 ; fi 2 )] 



2 » 



(3.11) 



constitutes a renormalization group invariant {i.e., //-independent) quantity [471 ] . In Fig. [7] 
we plot the combination above for the three different values of /i chosen; evidently the 
product of g 2 (fi 2 ), A(q 2 ;fi 2 ) and [1 + G(q 2 ; fi 2 )]~ 2 constructed from the solutions obtained 
is /i-independent to an extremely high degree of accuracy. 

The curves plotted for G(q 2 ) = u(q 2 ) (Fig. [6] left panel) should be compared with those 
obtained by Sternbeck [31] (reproduced in Fig. [8]), where the function (I2.2ip was studied 
in terms of Monte Carlo averages, and its asymptotic behavior was inferred from the iden- 



tity (12.201) : in both cases, the behavior is clearly the same. Note that in 31[ the extrapolation 
towards the zero limit was problematic, due to a lack of knowledge of the function L(q 2 ) 
[there denoted by q 2 v(q 2 )]; our analysis does not suffer from such a limitation, given that 
L(q 2 ) is completely determined by its own equation. 
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FIG. 8: The determination of the KO function obtained in [31]. Left panel: the function 
— u{p 2 ) = 1 — -F -1 (p 2 ), determined from Eq. (|2.20p by neglecting the L(p 2 ) component, plotted 
at three different renormalization points. Right panel: Comparison between the function —u(p 2 ) 
measured directly on the lattice and its asymptotic behavior —u(p 2 ). 



Now, the important point to emphasize is that the function G(q 2 ) = u(q 2 ) is not a fi- 
independent quantity; in fact, as we have established, its /i-dependence is exactly what is 
needed in order to enforce the //-independence of the RG- invariant expression given in ( 13. lip . 
In fact, its value at q 2 = 0, i.e. the KO parameter u(0), depends on the renormalization point. 
Notice that in the case of an IR divergent ghost dressing function the possible /i-dependence 
would be inconsequential, since, due to the identity (I2.20p . u(0) = —1 irrespectively of the 
value of fi chosen. Evidently, the situation is different in the case of an IR finite ghost 
dressing function, since u(0) acquires a non-trivial dependence on the renormalization scale. 
This dependence is plotted on the left panel of Fig. [9] for values of fi varying between 
2.6 and 4.3 GeV (due to the limited number of UV lattice points of our data), where we 
see that —G(0) = —u(0) varies in the interval [0.65, 0.68], in good agreement with the 
prediction (I2.34p . The dependence of w(0) on /i appears to be moderate, probably due 
to the rather narrow region of allowed [i values considered. A more detailed study with 
data sets extending deeper in the UV should allow one to explore the full extent of this 
dependence. 

On the right panel Fig. [9] we plot finally the //-dependence of the horizon function (12.331) . 
Notice that both dependencies can be fitted with a function that is characteristic of a phase 
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n[GeV] 



i[GeVl 



FIG. 9: Left panel: The dependence of the KO parameter u(0) = G(0) on the renormalization 
point u. Right panel: Same as in the previous panel but for the horizon function. In both cases the 
the continuous red line represents the fit given by Eq. (|3.13j) with a\ = 0.63, 61 = 3.57, C2 = 0.025 
and 02 = 28.61, 62 = 3.25, C2 = 0.05 for -G(0) and (/i(0)) 7= o respectively. 



transition, namely 



G(0) = -a^-bt) 

>7=0 



(3.12) 



as shown by the continuous red curves appearing in Fig. [9] (the values of the fitting parameter 
are also reported there). 

IV. CONCLUSIONS 



In this work we have presented an indirect determination of the KO function from recent 
lattice data on the behavior of the QCD gluon and ghost propagators Q in the Landau 
gauge. The results obtained are in very good agreement with the original study of the same 
quantity presented in [3JJ . 

Of particular interest is the observed dependence of the KO function, and particular of 
its infrared value u(0), on the renormalization point /i chosen to renormalize the theory. 
The /^-dependence of u(0), mild as it may seem at first sight, is definitely there, as one 
would expect, given that the KO function u(q 2 ) is not a renormalization-group invariant 
quantity, i.e., it is not intrinsically fx- independent. In fact, the observed \i- dependence is 
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a(n z ) = 0.330 and p = 3.6 GeV 
g(u 2 ) ; 0.295 and n ; 4.3 GeV 
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FIG. 10: The gluon propagator A(q 2 ) of the BFM, renormalized at three different points: 
fi = 3.0 GeV (black curve), \i = 3.6 GeV (red curve) and /i = 4.3 GeV (green curve). 

really sizeable when contrasted with the impressive absence of any /^-dependence displayed 
by a genuinely /^-independent quantity given in Eq. (13. lip which was computed using exactly 
the same sets of lattice data. We hope that the present work will contribute to the study 
of the possible effects that renormalization may have on the quantitative predictions of the 
KO formalism. 

Let us take a closer look at the background-quantum identity given in Eq. (I2.12p . which 
relates the conventional gluon propagator A with the gluon propagator A of the BFM. 
Eq. (12.121) assumes that the corresponding gauge-fixing parameters, namely the £ of the 
and the £q used in the BFM to gauge-fix the quantum fields appearing inside the loops, are 
equal (£ = £q). In the Landau gauge, £ = £q = 0, due to the central equality of Eq. (12.261) . 
we have that 

u{q 2 ) - 



- 1. 



Af 



(4.1) 



Interestingly enough, this simple formula expresses the KO function in terms of two gluon 
propagators calculated in the Landau gauge of two very distinct gauge-fixing schemes, with 
no direct reference to the ghost sector of the theory. This observation opens up the possibility 
of deducing the structure of the KO function using an entirely different, and completely novel, 



approach. Specifically, one may envisage a lattice simulation of A 



48| : then, u(q 2 ) may be 



obtained from (14.11) by simply forming the ratio of the two gluon propagators. Given that 



A(0) is found to be finite on the lattice [26|, |27|, it is clear that, in order for the standard 
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KO criterion to be satisfied (i.e., u(0) = — 1), A must diverge in the IR. Needless to say, we 
consider such a scenario highly unlikely. What is far more likely to happen, in our opinion, 
is to find a perfectly finite and well-behaved A, which in the deep IR will be about an order 
of magnitude larger than A(0), furnishing a value u(0) ~ —0.6, namely what we have found 
in our analysis. In fact, one may turn the argument around: combining the results of this 



article with the lattice data for A 
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271 ] . one may use (14. ip to predict the outcome of the 



lattice simulation for A; our prediction for the case of SU(3) is shown in Fig. [TUl 
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